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ABSTRACT 

We present high-resolution interferometric Sub-millimeter Array (SMA) imaging at 890/zm (~ 2" 
resolution) of two millimeter selected galaxies - MMJ100015+021549 and MMJ100047+021021 - 
discovered with the Max-Planck Millimeter Bolometer (MAMBO) on the IRAM 30 m telescope and 
also detected with Bolocam on the CSO, in the COSMOS field. The first source is significantly 
detected at the ~ Her level, while the second source is tentatively detected at the ~ 4<r level, leading 
to a positional accuracy of ~ 0.2 — 0.3". MM100015+021549 is identified with a faint radio and 
iGband source. MMJ100047+021021 shows no radio emission and is tentatively identified with a 
very faint iGband peak which lies at ~ 1.2" from a clumpy optical source. The submillimeter-to- 
radio flux ratio for MM100015+021549 yields a redshift of ~ 4.8, consistent with the redshift implied 
by the UV-to-submillimeter photometry, z ~ 3.0 — 5.0. We find evidence for warm dust in this 
source with an infrared luminosity in the range ~ 0.9 — 2.5 x 10 13 Lq, supporting the increasing 
evidence for a population luminous submillimeter galaxies at z > 3. Finally, the lack of photometric 
data for MMJ100047+021021 does not allow us to investigate its properties in detail, however its 
submillimeter-to-radio ratio implies z > 3.5. 
Subject headings: galaxies: evolution — galaxies: high-redshift — galaxies: starburst 



1. INTRODUCTION 

Submillimeter (submm) wavelength blank-field surveys 
discovered a population of heavily dust-obs cured star- 
burst galaxies at high redshifts (e.g. iSmail et al.l 119971 : 
iBarger et all 119981: Hughes et all 11998ft. These submm 
galaxies (SMGs) contribute a large fraction of the co- 
moving infrared (IR) l uminosity density at high-redshift 
(|Le Floc'h et al.l l2005f ). and their clustering properties 
indicate they could be progenitors of the most massive 
galaxies at z < 1 tiBlain et al.l 120041 : iViero et al.1 120091 : 
I Aravena et al.l[2010h . The number counts and dust pro- 
duction rates of these massive, high-redshift starburst 
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galaxies actu ally place tight con strains on galaxy forma- 
tion models (jBaugh et al.ll2005() . 

The identification of optical counterparts to the 
submm sources has mainly been based on the identifica- 
tion of radio counterparts, detected through deep Very 
Large Array (VLA) 20 cm imaging. Thi s technique ex- 
ploits the local far-IR/radio correlation (|Condonl 11992ft 
and statistical arguments to claim an association between 
a closeby radio source and the submm source, localizing 
the SMG with an accuracy of ~ 1". The direct search 
for counterparts in optical wavebands is not practical due 
to their common faintness and the high number density 
of optical sources within the positional uncertainty. For 
about 70 — 80% of all submm sources precise positions 
can be deter mined from an ide ntification of faint radio 
counterparts (jlvison et al.H2007J) . 

However, the depth of current radio observations limits 
the identification of submm sources to z < 3, since due to 
the strong K-correction, the radio flu x (unlike the s ubmm 
flux) drops rapidly with redshift (Cari lli fc Yunl ll999). 
The only unambiguous way to localize the highest- 
redshift SMGs in the optical/near-IR is via submm in- 
tcrfcrometry. 

Radio-selected SMGs h ave been found to hav e a me- 
dian redshift of z = 2.3 ({Chapman et al.l I2005T ). How- 
ever, recent interferometric submm imaging of a flux- 
limited sample of SMGs has provided substantial evi- 
dence that up to ~ 30% of this population is likely lo- 
cated at z > 3 (jlono et al.l 120061 : lYounger et all 120071 
2009). These radio-dim SMGs appear to be inconspic- 
uous at optical wavelengths, likely obscured by large 
amounts of dust, and their nature can only be re- 
vealed in the IR jlWang et al.l 120071 iDannerbauer et all 
120081 iCowie et all 12009ft . To date, only five such high- 
redshift SMGs have been confirmed spectroscopically 
(jCapak et al.l 120081 : ISchinnerer et alll2008t ICoppin et all 
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Multi-wavelength photometry 




Band 


X b 


MM1 


MM14S 


MM14N 


Units 


B+ 


0.45 


< 21 


< 21 


24 ±7 


njy 


i+ 


0.77 


56 ±16 


< 35 


116 ±11 


njy 


K S 


2.2 


768 ± 77 


222 ± 160 


340 ± 140 


njy 


IRAC 


3.6 


1.3 ±0.2 


< 0.2 


0.20 ± 0.06 


MJy 


IRAC 


1.5 


1.9 ±0.3 


< 0.3 


0.32 ±0.11 


MJy 


MIPS 


21 




< 0.06 


< 0.06 


mjy 


MIPS 


70 


< 2.1 


< 2.7 


< 2.7 


mjy 


MIPS 


160 


< 15 


< 15 


< 15 


mjy 


LABOCA 


870 


16.4 ±1.8 


< 10.5 




mjy 


SMA 


890 


16.8 ±1.5 


8.5 ±2.0 


<6 


mjy 


Bolocam 


1100 


6.9 ±1.9 


3.8 ±19 




mjy 


MAMBO 


1200 


6.2 ±0.9 


4.1 ±1.0 




mjy 


VLA 


20 


35 ±10 


< 30 


< 30 


A*Jy 


a Measurements at 3.6 and 4.5 i 


mi based on 


the deep Spi 


tzer 



warm mission exposures. b Wavelength in microns, except for the 
VLA where it is given in centimeters. c Photometry in the MIPS 
24 lira bands not possible due to blending with a bright source. 



l200llDaddi et al.ll2009rI at lKnudsen et al.ll2010l) . 

In this paper, we report accurate astrometry of two 
bright SMGs in the COSMOS field whose radio to 
submm flux ratio indicates that are likely at z > 3. 
We assume a standard cosmology with Ho = 71 km s _1 
Mpc"\ n A = 0.73 and D, M = 0.27. 

2. OBSERVATIONS 

2.1. The COSMOS field 

COSMOS is the largest deep survey carried out with 
the Hubble Space telescope (HST) covering ~ 2 deg 2 
in the sky. Extensive imaging of the COSMOS field 
has been performed from the X-rays to the radio wave- 
lengths. This includes: complete optical/near-IR cov- 
erage in 22 broad and intermediate bands with several 
ground based observatories including the Subaru tele- 
scope, the CFHT, the UKIRT and the KPNO; IR imag- 
ing with Spitzer, including new deep 3.6 and 4.5 /^m 
images obtained as part of its warm mission; and ra- 
dio imaging with the VLA at 20 cm. For details in 
the opti cal/IR imaging and cat alogs of the COSMOS 
field see iCapak et all ([2001 and lllbert et~aT1 (pOOl . A 
complete description o f the radio imaging is given in 
ISchinnerer et al 



ftiono 
(2007) 



2.2. Sample selection 

In the course of the MAMBO 1 .2 mm survey of the 
COSMOS field (|Bertoldi et al.ll2007l) . fifteen sources were 
detected with S/N > 4, five of which do not have a 
significant radio counterpart ( < 30 /iJy, 3cr). From 
these five radio-faint mm sources, we selected two that 
were also significantly detected at 1.1 mm with Bolocam: 
MMJ100016+021549 and MMJ100047+021018. Here- 
after, we refer to the se sources as MM1 a nd MM14, 
respectively, following iBertoldi et~ai] (|2007l ). The de- 
boosted flux densities for MM1 were Si. 2mm = 6.2 ± 0.9 
mjy and S Llmm = 5.9 ± 1.9 mJy in the MAMB O and 
Bolocam maps, respectively (Bcrtoldi et al. 2007). This 
source was recently detected with the Large Bolometer 
Camera (LABOCA) with S s70f , m = 16.4 ± 1.8 mjy (Al- 
brecht et al., in preparation). The MAMBO and Bolo- 
cam deboosted flux densities for MM14 were Si. 2mm = 
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Fig. 1. — Postage stamps centered at the position of the SMA 
detections. Images are 14" X 14" in size. The images to the left 
show the SMA 890 fim emission at significance levels of -4, -2, 2, 
4, 6, 8 and IOct (see text). Solid and dashed contours represent 
positive and negative fluxes, respectively. The red 2" diameter 
circle represents the position of the SMA detection. 

4.1 ± 1.0 mjy and S i.i mm = 3.6 ± 1.9 mjy, respectively 
(jBertoldi et al.ll2~007D . however no significant emission is 
seen in the LABOCA image down to 10.5 mjy (~3<r). 

2.3. SMA 890 [im imaging 

SMA observations of MM1 and MM14 were carried 
out in 2009 January 10 and 2009 May 07, respectively, 
under good weather conditions (T250GHZ < 0.1). The re- 
ceivers have two sidebands, each with 2 GHz bandwidth, 
which when averaged yield a 4 GHz effective bandwidth 
centered at 345 GHz (A ~ 890 /im). Seven antennas, 
arranged in the compact (COM) configuration^ with the 
MAMBO source positions as phase tracking centers. 

The data were calibrated with the MIR package 
(jScoville et al.lll993h specially adapted for SMA data, us- 
ing the strong continuum source 3C273 (S345GHZ ~ 4.9 
Jy) as passband calibrator and Ceres (S345GHZ ~ 4.2 Jy) 
for primary flux calibration. The flux scale is estimated 
to be accurate within 20%. The quasars J0854+201 
(-3.1 Jy; ~ 24.2° away) and J1058+015 (~ 1.7 Jy; 
— 14.5° away) were observed every ~ 25 min for gai n 
calibration. Following lYounger et all (|2008l I2009L 120101 1 , 
we also performed hourly scans of a dimmer, but signif- 
icantly closer test quasar J1008+063 (-0.11 Jy; - 4.5° 
away) to empirically verify the phase transfer and esti- 
mate the positional uncertainty. The visibility data for 
MM1 showed good phase stability, however, about half of 
the data for MM14 had to be flagged due to bad phases. 

The calibrated visibility data were imaged using the 
AIPS software. We used the AIPS task IMAGR, which 
uses the CLEAN algorithm, and natural weighting to 
deconvolve the images down to \a in a box centered 
on our targets. This led to beam sizes of 2.55" x 1.86" 
and 1.96" x 1.86" and noise levels of 1.5 mjy and 1.95 
mJy beam -1 for MM1 and MM14, respectively. Finally, 
fluxes were measured with Gaussian fits using the JMFIT 
task, included in the AIPS software. 

3. RESULTS AND ANALYSIS 

3.1. SMA detections 

Fitting two-dimensional Gaussians to the SMA images, 
leads to flux densities of 16.8±1.5 mjy and 8.5±2.0 mJy 

9 half-power beam width ( HPBW) « 2"; see 

http://smal.sma.hawaii.edu/specs.html for more information. 
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Fig. 2. — A'-band 8" X 8" postage stamp images centered at the 
position of the SMA detections of MM1 (left) and MM14 (right), 
with _R"-band white contours overlaid at 2, 3, 4 and 5a levels. Green 
contours represent the HST F814W image smoothed with a Gaus- 
sian kernel of 3 pixels (0.15") at 2, 3, 4 and 5a levels. The red 
contours represent the 890 fira emission with levels as in Fig. [l] 
The X-band image for MM14 has also been smoothed with a Gaus- 
sian kernel of 3 pixels (0.45") to enhance the significance level. 

for MM1 and MM14, respectively. The MM14 detection 
is tentative since, as we shall see below, it could not 
be reliably identified with any significant counterpart at 
other wavelengths. 

The Gaussian fit indicates that both sources are un- 
resolved with a maximum deconvolved FWHM size of 
1.8" x 1.0" and 2.0" x 0.6", respectively, consistent with 
the sizes found for high-redshift SMG s, which are typ- 
ically unresolv ed at ~ 2" resolut ion dlono et al.l 120061: 
Younger et ail l2007t iWang et al.l 120071; lYounger et al.l 
2009). This also yields from the flatness of the real vis- 
ibility amplitudes as a function of the projected base- 
line length. The measured SMA position for MM1 is 
a(J2000) = 10 h 00 m 15.612 s , <5(J2000) = +02°15'49.00", 
with a positional error in the Gaussian fit of 0.09", while 
for MM14 it is a(J2000) = 10 h 00 m 47.329 s , 5(J2000) = 
+02 o 10'21.44", with a positional error in the fit of 0.15". 
This positional accuracy in the fit is consistent with the 
one expected for the beam and S/N of the observations: 
0.09"and 0.2"for MM1 and MM14 respectively From 
the comparison of the reference ([Browne et al.l 119981) 
and measured (this work) positions of the test quasar 
J1008+063, we find a positional uncertainty of 0.18". 
This, added in quadrature to the positional error in the 
Gaussian fit to the MM1 and MM14 images, gives a po- 
sitional uncertainty of 0.2" and 0.27", respectively. 

3.2. Multi-wavelength counterparts 

MM1. The SMA peak position coincides with the 
position of a ~ 3.5er radio peak (Fig. [T|). A bright 
and elongated sou rce with photometric redshift z = 1.4 
(jllbert et al.l 120091) is located at ^2.1 "north- west from 
the SMA position. The radio peak lies within 0.3" from 
the SMA position, however given the beam of the radio 
image (~ 2"), we do not discard that part of the emission 
comes from this bright optical source. From Fig. [5] (left), 
we identify a very faint if-band source, located at sw 0.3" 
from the SMA position, as the likely counterpart. The 
bright optical source strongly contaminates the Spitzer 
images, making it difficult to reliably measure the faint 
emission of MM1. We extracted photometry in the IRAC 
bands by subtracting this bright source based on the K- 
band image convolved with the IRAC PSF (Table [TJ. We 
do not attempt to extract photometry of this source in 
the 24 [im images. No emission is detected at 70 and 160 



[im. 

The deep if -band images show two peaks, separated 
by about 0.6", or a physical scale of ~ 4.3 kpc at z ~ 3 
(Fig. [2]). The fainter peak appears to be the one asso- 
ciated with the submm emission, suggesting a possible 
double system, simila r to the case of the SMG AzTECll 
([Younger et alj|2009f) . 

MM14- The radio maps do not show any peak close to 
the position of the SMA source down to a 3c level of 30 
/ijy. At ~ 1.2" to the north of the SMA position we find 
a faint optical source that appears diffuse and faint in the 
if- band (Fig. [H). This source has a li kely photometric 
redshift of ~ 3.4 ([Mobasher et al.ll2007ft . 

From Fig. [I] the northern optical source appears to be 
composed by several "clumps" that extend over ~ 1.5", 
or ~ 11 kpc at z ~ 3.5. We also find a very faint K- 
band emission peak (^3cr in the smoothed image), lo- 
cated 0.4" to the south of the SMA position (Fig. EJ. 
This peak has ~ 2<r significance in the original if -band 
image (without smoothing; Fig. [IJ, implying the source 
is spatially extended. Due to its proximity to the SMA 
position, it appears to be the most likely counterpart 
to the submm emission despite its faintness. Hereafter, 
we refer to this source as MM14S (south), while for the 
northern optical source we refer as MM14N. However, 
since we could not reliably identify any significant multi- 
wavelength counterpart for this source and given the rel- 
atively low significance of the SMA detection, we label 
this source (MM14S) as a tentative detection. 

The spatial configuration between MM14S and 
MM14N could resemble a merger/interaction system, 
where the submm emission comes from a highly obscured 
source (MM14S), but could also correspond to an ex- 
tended galaxy with the submm emission located in an 
obscured spiral arm. The offset between the MM14N 
and the submm position (MM14S) is ~ 9 kpc, assuming 
z ~ 3. 5, similar to the case of the high-redshift SMG 
GN20 ([Iono et al.l|2006f ), where the submm and the op- 
tical peaks are separated by ~ 0.8", or ~ 6 kpc. Based 
on the local density of sources with z > 3, n — 0.002 
arcsec -2 , we find that the probability that a z > 3 op- 
tical source brighter than MM14N is located by chance 
within a distance of 1.2" from the SMA position is only 
P = 0.9%, thus supporting a physical association. 

O ptical/IR photometry w as performed with SExtrac- 
tor (|Bertin &; Arnoutsl [i996') in a 2" aperture, using the 
if -band images for detection. None of our targets was de- 
tected in the Spitzer MIPS bands, and we thus provide 
3(7 upper limits based on their local noise level within 
one beam. The measured flux densities at several wave- 
lengths for MM1, MM14S and MM14N are listed in Table 

3.3. Photometric redshifts 

MM1 a nd MM14S were not detected in the COSMOS 
catalogs ([Capak et alj|2007t lllbert et all 12009ft and thus 
there is no previous estimate for their redshift. 

Assuming that the correlation between the far-IR and 
radio integrated emission from local starburst galaxies 
also holds for distant galaxies, the subm m-to-radio flux 
ratio can be used as a measure of redshift (|Carilli fc Yunl 
11999ft . Figure |31 shows the 890 /im to 20 cm flux ratio as a 
function of redshift for the prototypical starburst galaxy 
Arp220, and illustrates the uncertainties introduced by 
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Fig. 3. — Redshift estimation based on the submm-to-radio flux 
ratio. The solid line represents this ratio computed for the IR lu- 
minous starburst galaxy Arp220 (which roughly has 7d us t = 45 
K). The filled area shows the flux ratio obtained for similar mod- 
els with dust temperatures from 30 to 60 K. Dashed horizontal 
lines show the flux ratios measured for the radio-identified submm 
sources with SMA detections (Younger ct al. 2007. 2009), and for 
thehigh- redshift SMGs HDF 850.1 and GN20 llCowie et al.ll2009l : 
Iono ct al. 2006). The upper pointing arrow represents the lower 
limit in the flux ratio for MM14S. 



different dust temperatures. From this, we derive a red- 
shift of 4.8+14 for MM1, and a lower limit z > 3.5 for 
MM14S. 

Given the lack of optical/IR photometric information 
for MM14S, we do not attempt to derive the redshift and 
dust properties of this source, and hereafter we focus on 
MM1 and MM14N. 

We derived optical/IR photometric r edshifts for MM1 
and MM14N using the HyperZ code (jBolzonella et al.l 
2000). We allo wed for a redshift range z = — 7, and 
using a lCalzetti et al.l (|2000T ) extinction law with an ex- 
tinction range Ay = — 4. For comparison, we used 
two different set of templates librar ies: 5 solar metal- 
licity stellar population models from lBruzual fc Charlotl 
(2003), similar to the Hy perZ defaults; and 10 templates 
from the SWIRE library (jPolletta et al.ll2007f ). including 
models for an elliptical, three spiral, four starburst and 
three AGN dominated galaxies. 

As shown in Fig. [4j the x 2 distribution and redshift 
solutions for MM1 and MM14N obta ined using both set 
of tem plates are similar. Using the Br uzual fc Charlotl 
(2003) templates, the best fit for MM1 is produced at 



= 3.1 



+0.5 



by a spiral galaxy model, with Ay = 1.4 
(black curve); while the best fit for MM14N is produced 
at z = 3.7+Q3 by a single burst model, with Ay = 1.0. 
Using the SWIRE library, the best fits are produced by 



a late spiral galaxy at z = 3.0_Qg and a QSO template 
at z = 3.6t° 2 3 for MM1 and MM14N, respectively. The 
quoted uncertainties correspond to the 90% confidence 
level derived from the fitting rou t ine. 

Each of the lBruzual fc Charlotl (|2003l ) templates follow 
a prescription for the stellar mass to luminosity ratio (e.g. 
a s a function of a ge). Using the best fits and assuming 
a iChabrieil (|2003l) initial mass function, we find stellar 
masses of ~ 1.0 x 10 11 M Q and ~ 3.6 x 10 10 M , for 
MM1 and MM14N, respectively. We note, however, that 
these stellar mass estimates are rough given the large 
uncertainty in redshift. 
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Fig. 4. — Photometric redshift results. Top: the x 2 values as 
a function of redshift obtained by fitting the optical photometry 
with the the Bruzual & Chariot (2003) templates (black curves) 
and the SWIRE templates (red curves) for MM1 (left) and MM14N 
(right); and by fitting the UV-to-mm photometry with the SWIRE 
templates (green curves) for MM1. Bottom: SED of MM1. The 
colors represent the different template libraries used. The gray 
curve shows the best fit dust model to the far-IR data at z = 3 — 5. 



For MM1, we also computed redshifts by fitting tem- 
plates to the UV-to-mm photometry. We used the same 
SWIRE templates, redshift and extinction ranges. The 
best fit is produced by an Arp220 model, with Ay = 0.3 
and z — 4.5. The value of x 2 is dominated by the ratio 
between the submm to optical light, indicating a possible 
range of redshift from ^3.0 to 5.0. Our redshift estimates 
are in good agreement with the redshift implied by the 
submm to radio flux ratio and supports a higher redshift 
than that obtained purely from the optical photometry. 

3.4. Dust properties 

As illustrated in Fig. H the far-IR SED of MM1 can be 
described by a modified black body spectrum. We fitted 
a model of the form S v oc -Bj,(Xdust)(l — e _T ")i where 
BvlTdust) is the Planck function. The opacity, t„, is 
proportional to the dust mass M^ust an d to v? , with the 
cmissivity index /3 = 1.5. Given the large uncertainties 
in redshift, we consider a range of redshift z = 3 — 5 to 
model the photometric data (70 — 1200/im). 

We find dust temperatures and masses in the range 
T dust = 50-65 K and M dust = 1.2-2.6 x 10 9 M s . These 
models imply far-IR luminosities Ls-iooOjum = 0.9 — 2.5 x 
10 13 L(7). Using the computed IR luminosities and the 
iKennicufrq (|1998l ) empirical calibration, we compute star- 
formation rates (SFRs) of 1800 - 5000 M Q yr" 1 for this 



4. DISCUSSION 

Our detection of two out of the five significant radio- 
faint sources discovered in the MAMBO survey of the 
COSMOS field supports the increasing evidence for a 
population of luminous SMGs at z > 3. It has recently 
been suggested that these galaxies may be the progeni- 
tors o f quiescent, massive galaxies at z ~ 2 f Capa k et al.l 
120081; iCoppin et al.1 ,2009. 2010|) . The old stellar popula- 
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tions found in these z ~ 2 quiescent galaxies require to 
have formed at z > 3 in short bursts of star formation. 

If we assume that four of these five MAMBO sources 
with no radio significant radio counterparts are real and 
located at z > 3 (e.g., we expect that one of the five is 
a fake detection due to flux boosting), it implies that up 
to ~ 30% (4/15) of the S250GHZ > 4 mJy detections are 
located at z > 3. This is well justified by the confirma- 
tion of MM1 and tentative detection of MM14 with the 
SMA, and of two other MAMBO sources with LABOCA 
(Albrecht et al., in preparation). Assuming that 2 — 4 
of the MAMBO sources are uniformly distributed over 
the comoving volume spanned by z = 3 — 5 within the 
area of the MAMBO survey, we find a surface density of 
^13 — 30 deg~ 2 and a volume density of ~ 1 — 6 x 10 -6 



Mpc 3 . This is consistent with the upper limit of > 10 6 
Mpc -3 (> 6 deg~ 2 ) from the spe ctroscopically confirmed 
z > 4 SMGs (jCoppin et al.ll2009T > and with the density of 
z ~ 3 — 5 systems with baryonic masses > 10 11 M& f rom 
predictions of theoretical models ([Baugh et a l. 2005). 
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